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Abstract 

We propose quite a new method of analyzing the dynamical chiral symmetry breaking in gauge 

theories. Starting with the non-perturbative renormalization group equation for the Wilsonian 

O' 

^SJ ■ fermion potential, we define the weak solution of it in order to mathematically authorize solutions 
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pi |. with singularity. The weak solution is obtained uniquely and it successfully predicts the physically 
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correct vacuum, chiral condensates, dynamical mass, through its auto-convexizing power for the 
effective potential. Thus it works perfectly even for the first order phase transition in the finite 
density QCD. 
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I. INTRODUCTION 

Among various methods to analyze the dynamical chiral symmetry breaking, the non- 
perturbative renormalization group is quite effective since it may include non-ladder di- 
agrams which cure the gauge invariance problemja]. In the lowest order approximation, 
the Wilsonian effective action is expressed by a scale dependent fermion potential V^^{x,t), 
where x = j d'^xil){x)il){x) is an operator variable, t = log(Ao/A(t)) is the renormalization 
scale. This potential keeps all information of the fermionic interactions, the free energy Vq, 
the dynamical mass Md, the 4-fermion interactions G etc as follows: 
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Vw(x, t) = Voit) + Mrf(t)x + Git)^ + ■■■ . (1) 

The renormalization group equation for G (in case M = 0) is written as 

which makes G diverge at a finite scale t^ for upper critical initial value. This blowup nature 
itself is a correct behavior since it expresses the divergence of the susceptibility due to the 



spontaneous chiral symmetry breakdown[ 



QQ. 



However, due to this divergence we can not 



go beyond tc, and there is no way to calculate infrared physical quantities such as the chiral 
condensate or the dynamical mass. 

II. PARTIAL DIFFERENTIAL EQUATION AND ITS WEAK SOLUTION 

Various methods have been used to bypass this singularity, e.g., the bare massj^, auxiliary 
fields [l|-l3|, etc. Here we propose a new direct method of solving the renormalization group 
equation as a partial differential equation (PDE)J7[: 

OM^ . »M . 0. M ^ ^. KM. ^ -^e-« mil . M'V']. (3) 

where M{x, t) is called the mass function. The 4-fermion interaction G corresponds to ^ 
at the origin, and its divergence at finite tc means that the above PDE has no global classical 
solution beyond tc- 

The expected behavior of physically meaningful solution M{x,t) is plotted in Fig. 1(d), 
where M has a jump at the origin after tc. In order to accommodate such singular solution 
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FIG. 1: NJL plot, (a) Characteristics, (b) Mass function, (c) Discontinuity, (d) Weak solution. 



with discontinuity, we write down the weak version of the PDEp, |7| 
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dx 
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dx M(x, 0)(/?(x, 0) = 0. 



(4) 



The weak solution ^ is defined as to satisfy the above equation for any smooth and bounded 
test function ip{x,t). The weak solution satisfy original PDE except for discontinuity and 
position of discontinuity x = S{t) is controlled by the Rankine-Hugoniot (RH) condition, 

dS{t) 



dt 



-[M+-M_] = /(M+)-/(M_), 



(5) 



where M_|_,M_ are right and left limit at the discontinuity point respectively. 

To obtain the weak solution, first we set up characteristic curves x = X{t) representing 



^ The authors greatly appreciate helpful comments by Prof. Akitaka Matsumura who told us how to 
construct the weak solution. 



the contour lines of M, which satisfies 

dX{t) df{M,t) 



(6) 



dt DM ■ 

Where these curves are crossing each other, we must pick up one curve and introduce discon- 
tinuity according to the RH condition, and then we get the unique function M{x,t). Fig. 1 
shows the Nambu-Jona-Lassinio (NJL) model example {g = 1.005(7c) of these procedures. 
The contour lines are plotted in (a), where at some finite tc they start crossing with each 
other. The cross section of the total contours is seen in (b), where the derivative ^ diverges 
at the origin at some tc- 

Note that in obtaining M{x, t) as in (b) there is nothing singular, and it is just a motion 
of 'string'. However, to get the renormalized potential we have to define M{x, t) as a unique 
function of x. Then the RH condition determines the discontinuity as shown in (c), and the 
mass function in (d). This is the weak solution of the PDE and it defines global t solution 
uniquely. 

III. WEAK SOLUTION RESULTS FOR THE PHYSICAL QUANTITIES 

We show results in the finite density NJL where the first order phase transition occurs. 
Contours and mass function are plotted in Fig. 2, where at the central region five-fold 
structure appears corresponding to the three-fold local minima. In the renormalization 
procedure, two discontinuities appear pairwisely, move towards the origin, and finally merge 
into one at the origin as shown in Fig. 2(a). 

Snapshots in the course of renormalization are shown in Fig. 3, where the mass function 
M{t,x), the Wilsonian fermion potential V\v and the Legendre effective potential for a 
are plotted. It is astonishing that our method of weak solution uniquely determines their 
singularity structures and the resultant Legendre effective potential is always convexized. 
This means the dynamical mass and the chiral condensates are uniquely calculated, and 
perfectly correct in the sense that even in case there are multi local minima, the lowest free 
energy minimum is always chosen automatically. This feature is quite a new finding and 
shows powerfulness of the purely fermionic non-perturbative renormalization group and its 
wealc so.ut,o„R. Th,s analysis has been apphed to QCD. even wth fin.te density or ,gn- 
ladder, and proved to work perfectly to give physical quantities without any ambiguityp. 
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FIG. 2: NJL g = 1.7gc, /^ = 0.7. (a) Characteristics and jumps, (b) Evolution of mass function. 
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FIG. 3: Evolution of physical quantities by weak solution (NJL g = 1.7gc, /U = 0.7, t = 0.01, 0.5, 
0.6, oo). (a) Mass function, (b) Wilsonian fermion potential, (c) Legendre effective potential. 
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